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We develop a semi-analytical approach to calculate the polarizability tensors of an arbitrary individual scattered The approach is 
based on the calculation of induced electric and/or magnetic dipole moments on the scatterer. By taking the advantages of the present 
approach, we calculate the individual polarizability tensors of an isolated scatterer in a homogeneous isotropic medium. Moreover, 
we obtain the polarizability tensors of scatterers located between two different isotropic media. Furthermore, due to the nature of 
the proposed method, we may determine the effective polarizability of a scatterer in a periodic array. To this end and for comparison 
reasons, we investigate two other methods for calculating effective polarizability tensors of scatterers in a two-dimensional periodic 
array. The proposed approach, in comparison to other reported approaches, is simpler, easily implemented, and does not require 
spherical harmonic expansion or complicated far-field calculations. We examine the validity of the proposed approach using several 
examples and compare the results with other approaches. 


Index Terms —polarizability tensors, multipole moments, effective polarizability, substrated scatterers. 


I. Introduction 

HE electromagnetic response of composite materials, in¬ 
cluding metamaterials, mainly depends on two important 
factors: 1) the properties of their “constituent elements” 0 and 
2) the manner in which these elements are arranged together. 
The latter determines the interaction effects between the ele¬ 
ments while the former demonstrates the general response of 
the composite. Indeed, the electromagnetic behaviour of the 
constituent elements may be described by the polarizability 
tensors 0. Polarizability tensors determine the relation be¬ 
tween the induced dipole moments and the external electro¬ 
magnetic field excitation for a particle with linear response to 
the fields 0. We may sufficiently explain the electromagnetic 
behaviour of many small particles by their dipolar moments 
or equivalently their polarizability tensors. Therefore, various 
analytical and numerical approaches have been employed to 
calculate the polarizability tensors of different particles l3l- 
0 - 

For instance, the polarizability of a dielectric sphere is 
analytically calculated in 0. Moreover, analytic formulas 
are derived for polarizabilities of wire dipole 0, chiral 0, 
and U-shaped particles 0. Furthermore, the polarizabilities 
of two nonreciprocal magnetic particles; i.e., Tellegen-omega 
and “moving”-chiral are formulated in 0. The proposed 
analytical approaches are not only difficult but also limited to 
simple structures. Moreover, small modifications to the particle 
structure result in a different analytical model. 

Consequently, for more complicated shapes, different semi- 
analytical and numerical methods have been proposed Q- 
q, uni. For example, Alaee et al. derived the polarizabil- 
ties of a planar omega structure using multipole moments 
approach 0. Their method was based on the expansion 
of spherical harmonics of the Mie coefficients lUTll in the 

'Also, called “inclusion”, “scatterer”, or “particle” in the current 
manuscript. It is called “meta-atom” in metamaterials. 


Cartesian coordinate system. One of the disadvantages of 
this method is the cumbersome calculations of the spherical 
harmonics. Furthermore, the method is only valid for particles 
in homogeneous host medium d. 

Recently, another numerical method is developed to re¬ 
trieve the polarizability tensors of general bianisotropic par¬ 
ticles nm It directly calculates the induced dipole moments 
of an excited scatterer and does not require complicated 
calculations and harmonic expansions. In this method, one 
probes the scattered far-fields of an excited particle and relates 
them to the radiation effects of a pair of electric and magnetic 
dipole moments. This method is based on the assumption 
of predominantly dipolar resonances and neglecting higher 
order multipoles. Although this method is comprehensive, it 
is not suitable for particles in inhomogeneous host media 
(such as substrated particles). Moreover, with this method, it 
is impossible to directly calculate the effective polarizability 
tensors of particles in periodic arrays. In a periodic array, the 
effective polarizability tensors are parameters which take the 
effect of the interaction of particles into account. Indeed, using 
the proposed method in lH3l one needs to separately calculate 
the interaction between the particles and then analytically 
connect them to the calculated individual polarizability tensors 
in order to determine the effective polarizability tensors. 

To overcome all the above deficiencies, we propose a 
method to calculate the polarizability tensors of a general 
bianisotropic particle. This method is based on the calculation 
of the near field response of a particle to the plane wave 
illumination. We use the induced local charges/currents on the 
particle to determine multipole moments and consequently the 
polarizability tensors. This way we provide an approach in 
which it does not only add the strength of different reported 
semi-analytical approaches in 0, lfl3fl 2 l but also enables one to 

2 That is, calculation of higher order multipoles as well as the polarizability 
retrieval of general bianisotropic particles. 
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determine the polarizabilities of particles located on substrates. 
Furthermore, with our approach we demonstrate how one may 
consider the case in which a particle is surrounded by similar 
particles. That is, how one can directly calculate the effective 
polarizability tensors of particles in planar periodic arrays. 

Hence, the paper is organized as following. The second 
section introduces our approach: the general ideas, formu¬ 
lation and numerical test set-up. It provides the procedure 
of evaluating polarizability tensors for an arbitrary isolated 
particle in a homogenous host medium such as free space. 
It is followed up by the considerations needed for substarted 
particles. Then, our approach is investigated to calculate the 
effective polarizabilities of particles in planar periodic arrays. 
Next, the concept of rotation matrix is introduced. At the end 
of this section, two more distinct methods for calculating the 
effective polarizabilities in arrays are presented for the sake of 
comparison. After presenting the analytical framework, in the 
third section, the proposed method is verified using various 
examples. 

In the third section, we start with the calculation of po¬ 
larizability tensors of an individual chiral particle in free 
space and then compare the results with the analytical method 
presented in a as well as the scattered far-held method 
presented in ED- Then, we calculate the so called effec¬ 
tive polarizabilities of these chiral particles when they are 
periodically positioned in a planar surface. Notice, the cal¬ 
culation of effective polarizabilities is impossible using the 
presented approaches in 0, ed since they use the scattered 
far-helds as their input data. We also compare the results 
with two other proposed methods introduced in section III-CI 
Next, we examine a non-reciprocal Tellegen-omega particle 
and compare its polarizabilities with those presented in lfl3l . 
Thereafter, we calculate the polarizability tensors of a nano¬ 
gold split ring resonator to demonstrate the capability of our 
proposed method in the terahertz frequency region. We also 
show that the higher order multipoles of such structure are 
negligible compared to its dipole moments. This is impossible 
to demonstrate with the presented method in fT3l . As the 
last example, we calculate the polarizabilitis of a truncated 
sphere on a substrate. This is impossible to calculate using 
the proposed semi-analytical methods in 0 , H 3 ) . The results 
are compared with the analytical methods presented in H4l . 
03 . We show that the presented analytical methods can 
only provide the electric polarizabilities of this structure and 
do not present the bianisotropic terms; whereas it is proven 
in 03 , 021 that the presence of the substrate may induce 
bianisotropy to the structure. Therefore, one needs to consider 
the bianisotropic terms in the calculation of polarizabilities as 
we perform in our approach. Finally, in the last section we 
conclude the study. 


II. Theory 

We first present a procedure for retrieval of the individual 
polarizability tensors of an arbitrary scatterer in a homoge¬ 
neous host medium. Then, by introducing the concept of 
rotation matrix, we extend the proposed procedure in order to 
be suitable for the substrated scatterers. Moreover, we expand 


our method in order to calculate the effective polarizability 
tensors of scatterers in planar periodic arrays. At the end of 
this section and aside the proposed method, we develop two 
different methods to calculate the effective polarizabilities of 
scatterers in planar periodic arrays. 


A. Extraction of the Polarizability Tensors 

The parameters which connect the induced electric and 
magnetic dipole moments p and m to the “local” electric 
and magnetic fields E loc and H loc are called polarizabilities. 
These relations may be summarized in the following compact 
forms 0: 


p 


"< 3 ee 

a Qm 


"e 1oc " 

m 


_< 3 me 

a mm 


H loc _ 


Here a ee , a mm , o; em , and a me are the electric, magnetic, 
magnetoelectric, and electromagnetic polarizability tensors, 
respectively. In a three-dimensional coordinate system, each 
of these polarizability tensors has nine components. Therefore, 
in order to completely characterize a particle response to an 
electromagnetic wave, there are 36 polarizability components 
that should be calculated. 

Let us first consider a particle in the origin of the Cartesian 
coordinate system in a homogeneous host medium with the 
intrinsic impedance 77 as shown in Fig. [l]-a. Now let us start 
with the calculation of a ^ 9 and a™ where i=x,y,z. According 
to Ref. fTOl and equation dTJ, we should determine the particle 
response (p^ and m^) to an x-directed exciting electric field. 
To create such an electric field, we apply two oppositely 
plane waves which are traveling in the z-direction with similar 
polarizations for their electric fields (green color vectors in 
Fig. [Q-a); i.e., E± = E 0 e j ^ tTkz \ where cc and k are the 
angular frequency and wavenumber, respectively [E Notice, 
± sign in the superscript refers to the propagation in in¬ 
direction. Since the two magnetic fields are out of phase, then 
summing these two plane waves copies a standing wave with 
an electric-field maximum and a magnetic-field zero at the 
place of the particle assuming the particle is small compared 
to the wavelength. As a result, the polarizabilities can be easily 
calculated as: 


^ee _ Pj + +Pj 
^ix 2E 0 


OL t = 


rry +m j 
l 2 E 0 1 


( 2 ) 


Here, the db sign in the superscript of the dipole moments, 
denotes for the induced moment by the corresponding plane 
wave in =bn-direction. In a similar way, subtracting these two 
plane waves creates a local maximum of the magnetic field 
with zero electric field. Therefore, one can write 


atT = 


Pi + -Pi 

2E 0 




= 


2E 0 


-77. 


( 3 ) 


We should mention again that the exciting field is assumed 
to be a constant value in the place of the particle. This is true 
when the size of the particle is sufficiently smaller than the 
operating wavelength. 

One may obtain 12 polarizabilty components from equations 
© and ©. Other components may be calculated similarly 
with E± = Eoe^^ and Ef = Eoett^W as the 


3 we consider the time dependence all over the manuscript 
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Fig. 1. Polarizability tensor extraction setup: a) an arbitrary scatterer is located in the origin of the coordinates and is illuminated by six orthogonal plane 
waves, b) an arbitrary scatter is located on a half space substrate, the six plane waves are rotated so that can excite the substrated particle, c) an array in the 
x-y plane that uses same setup as in (b). 


excitation electromagnetic wave (red color and brown color 
vectors in Fig. QJa). So far, we expressed that six orthogonal 
plane waves are needed to calculate all polarizability compo¬ 
nents. Moreover, we built the relations between polarizability 
components and the dipole moments. In the next step, we 
present how the required dipole moments may be calculated 
from the induced currents and/or charges lf20l . 

When an electromagnetic field illuminates a particle, the 
equilibrium of charges (either free charges or bounded ones) 
in/on the particle is changed. This excited particle can be 
modelled as an expansion of multipole moments. The first four 
successive electric multipole moments in terms of the induced 
charge density p(r) are expressed as ll20l 

q = fv p( r )dv, Pi = fy np(r)dv, 

(4) 

q,j = fy rirjp(r)dv. 

Where V is the volume of the particle and r is the position 
vector in the Cartesian coordinate system. Subscripts i,j,k 
denote components of a Cartesian tensor, and q, p i? and q- 
are the electric monopole, dipole, and quadrupole moments, 
and so forth, respectively. Since the total electric charge of the 
particle after the illumination remains zero, we have: q = 0. 
Similarly, the first two magnetic multipole moments are ex¬ 
pressed in terms of the induced current density J(r) as: 

mi = \ fv( T X J( r ))i rfv > 

m ij =§fv( rx J(r)hrjdv, 

where mi, and my are the magnetic dipole and quarupole mo¬ 
ments, respectively. Finally, we may summarize our procedure 
as following. First, we separately subject the scatterer to three 
couples of plane waves in three orthogonal directions. Then, 
using the dipole moments calculated from equations @ and 


(0) we may easily find the polarizability components from 
equations 0 and 0. 


B. Rotation Matrix 

With the proposed method we require the wave propagation 
in three orthogonal directions to determine the whole polariz¬ 
ability tensors of a particle. Therefore, this approach can be 
easily applied to any arbitrary scatterer in a homogeneous host 
medium. However, in the case of substrated scatterers (Fig. [T]- 
b) or scatterers in planar arrays (Fig.Ec), one cannot attribute 
a specific wave number to the waves propagating along the 
interface of the substrate or the array. Therefore, the interface 
plane is forbidden in our methodology and one may not simply 
use equations 0 and 0 for the waves propagating in the 
directions along the interface; that is, x and y in Figs. [T]- 
b, Ec. However, since the polarizabilities of particles are 
characteristic parameters of particles, then they should not 
depend on the propagation direction. As a result, we propose 
to rotate the excitation fields so that there is no incident wave 
along the interface. 

Now, let us first assume the particle or the array is located 
at the interface of two different media with the intrinsic 
impedances 77 + and as shown in Figs. 0 b or E c. The 
interface is assumed to be at the xy- plane. We then rotate the 
coordinate of the incident waves and find the polarizability 
tensors in the rotated coordinate system (x f — y' — z'). After 
that, we perform the same procedure for two traveling waves in 
the opposite directions =b z'\ i.e., E f ± = Eoe ^ utTkTz ‘^ (green 
vectors in Fig. Eb). This way we obtain the polarizability 
components as we did in equations 0 and 0 , but now in the 
rotated coordinate system; i.e., a jj, a™, a\™, and a™? with 
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( 6 ) 


Notice, equations © provide 12 polarizability components of 
a/an particle/array when located between two different media. 
The remaining components can be calculated in a similar way 
using E'± = Eoei^ tTk±x '^ and E'± = i? 0 e J 7 a;t= F /c= 'V) as 
the excitation fields (red and brown colors in Figs. |T]-b or [T]- 
c). The next step is to convert these polrizability components 
from x' — y' — z' to x — y — z coordinate system. Indeed, the 
rotation matrix which defines the relation between the rotated 
and original coordinate systems may be found as ifTTl 


(7) 


where 
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and, 0 and <f are the rotation angles around the x- and y- axes, 
respectively. Finally^ the polarizability tensors in the rotated 
coordinate system ( a ') can be transformed to the polarizability 
tensors in the original coordinate (<S) using the rotation matrix 
R by the following transformation (see in Appendix-lAt: 


a = R • a • R. 


( 8 ) 


Where R~ x is the inverse of the rotation matrix R. 

The next point is associated with the polarizability tensors 
of a particle in an array. Since we take the advantage of the 
local currents in our approach, then we inherently take the 
interaction effect of other particles in the array into account. 
Consequently, when we use our approach for a particle in 
an array, we directly calculate the “effective” (or collective) 
polarizability tensors |[22l of the array. However, in order to 
compare our method with other approaches, we present two 
alternative methods for calculation of the effective polarizabil¬ 
ity tensors. The first one is the calculation of effective polar¬ 
izability tensors using the individual polarizability tensors and 
the interaction constants. The second method is the retrieval 
of the effective polarizabilty tensors from the reflection and 
transmission coefficients of the array. These two methods are 
described in the following section while an example which 
shows the comparison between these approaches is given in 
section 


C. Effective Polarizabilty Tensors of Planar Periodic Arrays 

1) Effective Polarizabilities from the Individual Polariz¬ 
abilities and the Interaction Constants 

In order to analyse a scatterer in the close proximity of 
other scatterers, the effects of other scatterers may be replaced 
by the interaction constants. According to equation ©, the 
dipole moments are expressed in terms of the local electric 


and magnetic fields E loc and H loc . These are the fields acting 
on the particle. Therefore, in a periodic array, these local fields 
are the contributions of the incident fields E inc and H mc as well 
as the fields created by other dipoles. We may formalize these 
contributions as m 

E loc = E inc + pee . p 

H loc = H i„c + . m W 


where /3 ee and /? mm are the electric and magnetic interaction 
constants. For a periodic array, they may be expressed as |[2): 


/3 ee = 


AI 0 

o h 
0 0 


0 

0 




mm _ 


( 10 ) 


P\\ ~ “ jWR^) e ~ jkR ° ’ 

P± ~ -3% i( 1 + jFR^) e ~ jkR ° ■ 

Here a is the period of the array and Rq = a/ 1.438 l2l . 
Replacing © in (Q]), results in the dipole moments in terms 
of the incident fields; i.e., 


P 

* 

E inc " 

m 

= a • 

H ‘ nc _ 


( 11 ) 


Where a is the “effective” polarizabilty tensor and can be 
calculated as 

* _\I- <S ee • _/3 ee -S em • yS™ ' 

_ -s me • y e i- <5 mm • ^ mm 



a ee 

a em 


_<5 me 

a mm 


( 12 ) 


In summary, we first calculate the polarizability tensors of an 
“individual” particle and then converts them to the “effective” 
polarizability tensors through (IT2t using the interaction con¬ 
stants presented in (ITOl) . Note that the aforementioned analysis 
is true for those periodic arrays which can be modeled by 
dipole moments only. Moreover, it is assumed that the electric 
and magnetic dipole moments of all particles are in one plane 
and orthogonal so that the effects of the electric/magnetic 
dipoles on the magnetic/electric dipoles are zero; that is, 
/3 em = ^ me = 0. 

2 ) Effective Polarizabilities from the Reflection and Trans¬ 
mission Coefficients 

In this approach, we retrieve the effective polarizability 
tensors from the reflection and transmission coefficients of the 
array. For our purpose, it is enough to use the normal illumina¬ 
tion. However, the method can be similarly generalized to the 
oblique illumination. Let us first consider the array in the x — y 
plane as in Fig. [l]-c. Then, we illuminate the array with two 
electromagnetic plane waves with x-polarized electric fields 
from two opposite directions — z and +z. The reflected and 
transmitted electric fields can be expressed in terms of the 
dipole moments as lf23) 


E re f ,± = ± Z X m), 

E tmns,± = E inc,± _ j^( yp zp z X m), 


( 13 ) 
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where we distinguish between the ±z illumination direction III. RESULTS AND DISCUSSION 

by the ± sign. Substituting ffi) into ® results in A A Wjre Chlral Partide 


E L ef,± = T my) = 

+ wTyH™ T &™E™ T ^Z H 'T) = 

± &T y T - ia” m )4 nc , 

for the x-component of the reflected electric field. Therefore, 
the co-component of the reflection coefficient reads as 


One of the interesting bianisotropic particles with a variety 
of applications in the field of electromagnetics is the “chiral” 
particle lf24l - ll26l . A chiral is an object which can be distin¬ 
guished from its mirror image and cannot be superposed onto 
it (27) • A typical wire chiral particle is shown in Fig. [2]-a. It 
is composed of a planar loop connected to two straight wires 
orthogonal to the loop plane. According to the coordinates of 


iC 


Ef 

E L nc 




ee dz (T em =F d me 

xx xy i u 


yx 


O- (!4) 


Similar formulas can be derived for the co-component of the 
transmission coefficient and cross-component of the reflection 
and transmission coefficients as 


T ± = 


^itrans 


= 1 - j^(v & xx ± & Ty ± & yx 
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yyy b 






E y 

7> 


iritrans 

rp± _ E y 
cr E inc 


± A em zb a me + -a mm ) 

J 2 a 2 \ I yx yy xx ' y xy )"> 

± a Qm qp a mQ - -a mm ) 

j 2 a 2 \ I yx yy ' xx y xy r 



and (fl5l) in order to retrieve the eight unknown “effective” 



polarizabilies. The results read as 

^xx — — 2), 

*™=J^(ZT co -ZR co -2\ 


Fig. 2. The scatterers under study: a) A wire chiral with radius of a and arm 
length of l, the radius of wire is ro, b) A Tellegen-Omega particle composed 
of two crossed PEC wires with total length of 21 on top of a ferrite sphere 
with radius of a, the ferrite sphere is biased along z-axis, c) A split ring 
resonator with an arm length of L and the thickness of t, d) A truncated gold 
sphere on a half space Saffir substrate. 
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Fig. [2]-a, this chiral particle provides the following ten major 
polarizability components 


Q/ = 7 
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= j£;( AT cr + A ^cr), d™ = -j£(AT cr - A R a ). 

06) 

Where SA co / cr = A+ /cr + A co/cr and AA co / cr = A+ /cr — A co/cr . 

Other “effective” tangential polarizability components can be 
retrieved using y polarized incident waves. Moreover, the ef¬ 
fective normal polarizability components (a ] kl ; where k and/or 
l =z and i,j= e and/or m) cannot be found using the normal 
illumination. In order to find for the normal polarizability 
components, one needs to extend this approach for the oblique 
incidence. However, it is not achievable without increasing the 
problem complexity. 

In the next section, we focus on different examples to 
validate our semi-analytical approach in calculation of the 
individual or effective polarizabilities. In the beginning, we 
calculate the individual polarizability tensors of a chiral parti¬ 
cle in microwave frequencies. Then we continue to calculate 
the effective polarizability tensors of a planar periodic array of 
these particles using the local currents and compare the results 
with two other presented methods in sections III-C1I and III-C2I 


a mm ^ a™™ XX, 

(17) 

<3 me ^ <£xx + <;xy. 

Other polarizability components are negligible or zero com¬ 
pared to these components. Moreover, from the reciprocity we 
have 

cJ ee = ((J ee ) T , a mm = (a mm ) T , a em = — (a me ) T . (18) 

In Ref. 0, a wire and loop model is used to analytically 
calculate the chiral polarizabilities in terms of the model 
parameters. We choose a wire chiral example with the radius 
a = 1.7 mm and an arm length of l = 1.1 mm which provides 
a resonant response around 10 GHz. The polarizabilities of 
this particle are calculated using three methods: local currents 
(our proposed method), the wire and loop model presented 
in (3), and the scattered far-field method proposed in (T3) . 
The results are depicted in Fig. [3j As can be seen from Fig. 0 
the results of our method coincide with those of scattered 
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Frequency (GHz) Frequency (GHz) 

“ Real: Wire & loop model “““'Imag: Wire & loop model 
— Real: Local currents — imag : Local currents 
Q Real: Far fields Q Imag : Far fields 

"IS 25 35 

Frequency (GHz) 



Fig. 3. The major polarizability components of the chiral element shown in Fig. [2| which are calculated by the proposed method and compared to the wire 
and loop model as well as the far-held method. 


far-field method ED- Moreover, the wire and loop model 
results in an accurate prediction of the resonance frequency 
while it does not precisely foretell about the strength of the 
polarizabilities at resonance. The difference is due to the 
fact that a chiral particle is not exactly a combination of a 
wire and a loop as considered in is. Indeed, they interact 
together which is neglected in that model. Notice, all results 
respect the reciprocity theorem through Eq. |TH1 Note also the 
different resonant frequencies of af z with other polarizability 
components. This is due to the smaller resonant length of the 


particle in the z direction. 

Next, having these individual polarizabilities in hand, we 
calculate the effective polarizabilities of a planar periodic array 
of such chiral elements. 

B. A Planar Periodic Array of Chiral Particles 

We build a planar periodic array composed of similar chiral 
particles of section IIII-AI in the xy- plane. The period of the 
array is assumed to be 4 mm in both x and y directions. For 
simplicity and comparison reasons, it is enough to consider 
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Frequency (GHz) Frequency (GHz) 


Frequency (GHz) 


— Real: Interaction constants 
—— Real: Local currents 
Q Real: Reflection/Transmission 


----- Imag: Interaction constants 

- Imag: Local currents 

n Imag: Reflection/Transmission 


Fig. 4. The effective polarizabilities of an array of chiral elements with the period of P = 4 mm. In the proposed method we used integration of local 
currents over a unit cell of the array while in the interaction constants method we used the individual polarizabilities of Fig. [3] and calculated the effective 
polarizabilities using the equation dl2b . In the retrieving method, we retrieved effective polarizabilities from the reflection and transmission coefficient of the 
array through equations ns. 


only the normal illumination. With this illumination, we may 
extract all major tangential polarizabilities of the chiral particle 
given in ([171) . However, the normal polarizability af z may 
not be extracted with normal illumination. Thus, nine major 
effective polarizabilities from ten important polarizabilities are 
calculated with three proposed methods: I) our method of 
direct calculation of the dipole moments through the integra¬ 
tion from the local currents over a unit cell of the array, or 
briefly the method of local currents , II) the proposed method 
which we convert the individual polarizabilities calculated in 
section- HTl-AI to effective polarizabilities using the interaction 
constants through equation (fl2b . or briefly the method of 


interaction constants , III) the proposed method of reflection 
and transmission coefficients using equations (fl6l) . or briefly 
the method of reflection/transmisssion. The results are plotted 
in Fig. |4l It is obvious from the results, that the three methods 
are in quite good agreements. However, our first proposed 
method ( local currents) is fairly easier to apply since we only 
need to directly measure the local currents in a unit cell. More¬ 
over, we do not require approximate assumptions ( interaction 
constants ) or complex calculations ( reflection/transmisssion ). 
We should mention that the method of local currents is 
independent of the complexity of the problem and we may 
calculate all 36 “effective” polarizability components without 













































































assuming any approximations or complex calculations for the 
oblique incidence. For instance, in the current example using 
our first proposed direct method, we may easily calculate 
the last major polarizability component af z . However, the 
other two methods require either approximations or complex 
calculations. 

Next, we examine the applicability of the method of local 
currents in calculation of the polarizabilities of non-reciprocal 
scatterers. 

C. A Non-reciprocal Tellegen-omega Particle 

Nonreciprocal particles and their applications are com¬ 
prehensively discussed in Q. These particles may be 
realized using a biased magnetic element (e.g. ferrites) 
and some conducting wires. Two famous examples are 
Tellegen-omega/“moving”-chiral particles which simultane¬ 
ously presents the properties of Tellegen and omega/“moving” 
and chiral bianisotropies, respectively. A Tellegen-omega par¬ 
ticle as shown in Fig. [2]-b may be composed of a ferrite sphere 
and two orthogonal wires. When an incident electric field 
excites one of the wire arms, the induced electric current will 
produce two crossed magnetic moments in the ferrite sphere 
which in turn they induce currents on the wires. From the 
symmetry of the structure with respect to the x and y axes, 
one can deduce that 

ee/mm/me/em _ ee/mm/me/em 

cx, xx ^yy ’ 


^ee/mm/me/em _ ^ee/mm/me/em 
.'Ey Uj yx 

We calculated the major polarizability components of the 
structure using the presented method of local currents and 
compared them with the method of scattered far-field presented 
in m. The results are shown in Fig. \5\ Again, for this 
non-reciprocal element, the results of our proposed method 
coincide with the scattered far-field results. This proves the 
power of our method in calculation of the polarizabilities of 
non-reciprocal particles. 

The next example is associated with the calculation of the 
polarizabilities of an artificial optical magnetic dipole; that is, 
a nano-gold split-ring resonator (SRR). 

D. A Split Ring Resonator 

Split ring resonators (SRRs) are widely used in terahertz 
frequencies since they provide artificial magnetism which is 
naturally prohibited at these high frequencies ll28l . A typical 
SRR is depicted in Fig. Oc. Its magnetic moment is created 
in y direction due to a circulating current in the SRR arms. 
This current may be induced by an electromagnetic plane wave 
illumination. 

The main polarizability components of the proposed SRR 
are calculated using our method of local currents and are 
compared with the method of scattered far-field presented 
in ca. The results are shown in Fig. [6] As can be seen from 
the figure, the results of the proposed method coincides with 
the results extracted from the approach of the scattered far- 
field. Therefore, one may deduce that the SRR can be replaced 
with an electric and a magnetic dipole. 


x 10" 19 
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Imag : Local currents Q Imag : Far fields 

Fig. 6. The Major polarizabilities of the proposed split ring resonator. 
Comparison of our method with the far-field method. The length and the 
thickness of the SRR are L = 300 nm and t = 40 nm while it is made of 
gold in which its dispersive material parameters are taken from (29J. 


At this step, we clarify an important issue concerning the 
calculations of higher order multipoles. It is impossible to in¬ 
vestigate the effect of higher order multipoles using the method 
proposed in Ref. fT3l from the scattered far fields. Moreover, 
to calculate higher order multipoles by the expansion of 
spherical harmonics using the scattered fields one requires 
to solve a complicated problem as the order increases (e.g. 
see in Ref. C2). However, since in our proposed method we 
are numerically integrating from the induced charges/currents 
over the particle, then we may directly calculate higher order 
multipoles up to any desired order without increasing the 
problem complexities. A common way to investigate the 
effect of higher order multipoles is to compare the effect of 
different multipole moments in the scattering cross section 
of the particle lfl2l . As the orders of multipoles increase 
their contributions in the scattering cross section decrease. 
Here, we compare the effect of strongest electric (p x ) and 
magnetic dipole moments (m y ) with the strongest electric 
quadropole moment (Q xz ) in Fig. [71 As can be seen from 
this figure, the effect of quaropole moments are negligible 
and therefore SRR can be effectively modelled as a pair of 
electric and magnetic dipole. This is an important fact that 
we may investigate directly using our proposed method of 
local currents. This way we may have a fast physical insight 
of the complexity of a desired particle. Indeed, without this 
analysis, it is impossible to consider a particle as dipoles 
only. As a result, our proposed method not only provides the 
polarizabilities of any kind of dipole particles but also present 
a very useful and powerful physical insight of any general 
particle. More importantly, we do not need extra calculations 
to grasp such physical insight. Therefore, we gain too many 
information with a simple method. 

Next, we present the potential of our proposed method for 
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Fig. 5. Major polarizabilities of the Tellegen-Omega particle of Fig. [2)b calculated through proposed method and compared to far-held approach. The 
properties and dimensions of the particle are: l = 1mm, a = 0.5mm, Ferrite material: YIG with relative permittivity e r = 15 and magnetic saturation 
M s = 1780 G. The applied bias held is 3000 Oe. 


particles supported by a refractive substrate. 



Fig. 7. Scattering cross sections due to the different multipole moments of 
the SRR. 


E. Particles on Substrates 

It is very famous, practical and essential to find the polar¬ 
izabilities of a substrated particle fTTl-HTbil. Indeed, in most 
practical cases the particles cannot be freestanding and should 
be supported by a substrate for mechanical robustness. The 
substrate effect on the polarizabilities of a particle is inves¬ 
tigated in many works G3-IH3- Indeed, it is impossible to 
apply the proposed scattered far-field method presented in lfT3ll 
to calculate for the polarizabilities of a substrated particle. 
There are some limited analytical works (see e.g. in H4l . 
fT~5l0 which solved this problem for a substrated gold trun¬ 
cated sphere. They have solved laplace equations and applied 
appropriate boundary conditions. Although in their works the 
electrical polarizabilities in the co and cross directions were 
calculated, three other polarizabilities (i.e., magnetic, magne¬ 
toelectric, and electromagnetic) were neglected. We show here 
that some of these components are comparable to their electric 
counterparts. This is due to the presence of the substrate 
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Fig. 8. The electrical polarizabilities of a gold truncated sphere on a Saffir 
substrate shown in Fig. |2jd with ^ = —0.4. 


which may lead to a kind of bianisotropy called substrate- 
induced bianisotropy (SIB) ifTbl . Therefore, we calculate all 
polarizabilities of a substrated particle to examine the effect 
of SIB. We apply the proposed method of local currents to 
calculate the polarizabilities while we consider the concept 
of rotation matrix introduced in sectioi fH-Bl to the structure 
presented in Refs. Ifl4l , JT5). This structure is shown in Fig. [21- 
d. It is a truncated gold nano-sphere on a Saffir substrate 
(e r = 3.13). We create different geometries by changing the 
truncating ratio (^). The negative and positive ratios are related 
to the island and cap shapes, respectively. When ^ = 0, we 
have a hemisphere and with ^ = —1, we form a complete 
sphere. As the structure is symmetric in the xy- plane, one can 
define normal and parallel polarizabilities as ay = a xx = a yy 
and a± = a zz . Moreover, all calculated polarizabilities are 
normalized to the volume of corresponding truncated sphere. 

The normalized parallel and normal electrical polarizabili¬ 
ties for ^ = —0.4 are plotted in Fig. [5] using three different 
method; that is, our local currents, the method proposed 
in rm and the presented method in (T5) . The results of 
these three methods show similar frequency behavior with 
some differences. The differences between the methods of 
Refs. fH and C3 are due to approximation used in their 
numerical integrations. For example, in M, it is stated that 
the results may be deviated up to 13.3 percent with changing 
the number of multipoles taken into account in particle- 
substrate interaction. However, the disagreements between our 
method and the other two methods is due to the fact that 
they have neglected the effect of bianisotropy induced by the 
substrate which we present in the following. 

As we discussed earlier, our proposed method can provide 
all polarizability components of the structure. With our anal¬ 
ysis, we observe that there is a magnetoelectric polarizability 
component in the xy plane which is comparable to the 
previously calculated electrical polarizabilities in Refs. H41 . 
na. Although it is small for small truncation ratio, it starts to 


10.0 



Wavelength (nm) 


Fig. 9. a) The magnetoelectric polarizability of a sphere with ^ = — 1 
compared to its normal and parallel polarizabilities, b) the same figure for a 
hemisphere ^ = 0. 


0.25 



0.051 ------ 
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Fig. 10. The ratio of magnetoelectric polarizability to the parallel electric 
polarizability for different truncation ratios. 


be significant as we increase the truncation ratio. For instance, 
it is about one tenth of the parallel electric polarizability 
component at the resonance for a complete sphere ^ = — 1 
(Fig. [9]- a) while it increases to one fifth of the parallel electric 
polarizability component for a hemisphere ^ = 0 (Fig. [9} 
b). Notice, the magnetic polarizability (not shown here) is 
much smaller than the plotted polarizabilities in Fig.[9]-a,b and 
may be correctly neglected. Figure. flQl demonstrates the effect 
of the truncation ratio on the magnetoelectric polarizability 
component. As a result, we prove that the substrate has 
certainly an inportant effect on the polarizabilitry components 
of a particle which should be considered in modeling. With 
this example we conclude our study on the calculation of 
polarizability components of different type of particles under 
different circumstances. 
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IV. Conclusion 
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We presented a method for the calculation of polarizabil¬ 
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with substrated particles using an example of truncated sphere 
on a substrate. Indeed, we investigated that the presence of 
the substrate lead to a substrate-induced bianisotropy which 
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Generally with our approach, we filled two incomplete parts 
in the field of calculation of the polarizability tensors; that is, 
firstly calculation of the polarizabilities for substrated particles 
and secondly calculation of the effective polarizabilities of a 
periodic planar array. 


Appendix A 

Transformation of polarizability tensors from 

ONE COORDINATE TO A ROTATED COORDINATE 


We introduced the rotation matrix in the manuscript through 
equation ©. Therefore, we can transform every vector in the 
rotated coordinate x'—y'—z' to the original coordinate x—y—z 
through A' = R(6, 0)A. Consequently, for the fields-moments 
relation in the rotated coordinate system; i.e., 


v 

in' 


= a 


E r 

H' 5 


one may write: 


R- 


P 

m 


a' -R- 


or equivalently: 


P 

m 


= R~ 1 - a' -R- 


E 

H ’ 


( 20 ) 


( 21 ) 


( 22 ) 


Note, vectors without prime sign are original vectors and with 
prime sign are vectors in the rotated coordinate system. Finally 
from (l22k one may extract: 

a = R- 1 ■ a' • R. (23) 
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as the relation between the polarizability tensors in the original 
and rotated coordinate systems. 














